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Abstract Granada oscillation code (GraCo) is a software 
constructed to compute adiabatic and non-adiabatic oscilla- 
tion eigenfunctions and eigenvalues. The adiabatic version 
gives the standard numerical resolution, and also the Richard- 
son extrapolation, different sets of eigenfunctions, different 
outer mechanical boundary conditions or different integra- 
tion variables. The non-adiabatic version can include the at- 
mosphere-pulsation interaction. The code has been used for 
intensive studies of d Scuti, 7 Doradus, j3 Ceph., SdO and, 
SdB stars. The non adiabatic observables "phase-lag" (the 
phase between the effective temperature variations and the 



radial displacement) and 



ST, 



'ff 



^'ff 



(relative surface temperature 



variation) can help to the modal identification. These quan- 
tities together with the energy balance ("growth rate") pro- 
vide useful additional information to the adiabatic resolution 
(eigenfrequencies and eigenfunctions). 
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1 Introduction 



GraCo (iMova et all l2003h is a software developed to solve 
non-radial adiabatic and non-adiabatic oscillation equations. 
It is written in fortran95 language. It can be used for models 
all over the HR diagram. GraCo is able to work with three 
differ ent sources of e quilibrium mo dels: CESAM (Morel ,_ 
1997h . Granada Code ( Claret Lfl999h and JMSTAR dLawlor am 
20061) . The numerical technique used is t he so called Henyey 
relaxation method as it is described in lUnno et alTI (Il989l) 
(Section 18.2). The simple representation of the system of 
differential equations in terms of second-order centred dif- 
ferences is adopted for the numerical resolution. 
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2 Adiabatic case 

T he adiabatic s ystern of differential equations is described 
in lUnno et al. I (Il989h (p. 161). The code has the possibility 
of choosing between two sets of eigenfunctions. Both sets 
include the radial displacement (^^) and the Eulerian pertur- 
bation of the gravitational potential (-(/)'). The sets differ in 
the use of the Lagrangian or the Eulerian variation of the 
pressure, and the addition or not of a function of the radial 
displacement {U^r/r, with U = dlnm/dlnr, m the mass and 
r the radius), to the derivative of the Eulerian perturbation o f 
the gravitational potential {d^' /dr) dVorontsof et al. Lll97q) . 

As b oundary conditions GraCo uses those prescribed in 
lUnno et al . (1989) (pp 162 j^. The solutions must satisfy 
regularity conditions at the innermost mesh-point. As a first 
surface condition the continuity of (/)' and its first derivative 
are imposed. As second surface condition, the mechanical 
one, the program offers two possibilities: I) The Lagrangian 
variation of the pressure vanishes {dp = 0), or 2) makes use 
of the isotherma l reflective wave boundary condition (see 
lUnnoetal. 1(11989 1. pp 163^^. 

Another degree of freedom of GraCo is the variable of 
integration. The program can solve the system of differential 
equations as a function of the logarithm of the radius (In r) 
or the ratio between the radius and the pressure (-p). The 
first one largely weights the inner regions and the second the 
outer ones. Depending on the physics to be tested, the user 
can choose the most convenient variable. 

In order to solve the eigenvalue problem, the outer bound- 
jare cendiiiop for the gravitational potential is removed. In 
this case, foi" each trial eigenfrequency we have a unique 
solution. But not all the solutions obtained with every trial 
eigenfrequency fulfill the removed boundary condition. The 
spectra of the system is the set of eigenfrequencies and eigen- 
functions fulfilling this outer boundary condition. 

As the system of differential equations has been replaced 
by a system of centred difference equations of second order, 
the truncation error in the eigenfrequencies and the eigen- 
functions are of the order N~^, with A^ the number of mesh 
points of the equilibrium model. To obtain more accurate 
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eigenfrequencies, the cod e can make use of the so called 
Richardson extrapolation (IShibahashi et al. L Il98lh . a com- 
bination of the values obtained with A'^ and N /2 mesh points 
cancelling the leading order error. 

Finally, GraCo can compute the firs t order rotationa l split- 
ting through the Ledoux coefficient dLedouxL Il95lh . The 
eigenfrequencies for radial modes can be obtained in two 
ways: 1) Using the LAWE second order differential equa- 
tion, or 2) setting ^ = in the standard non-radial system of 
equations. 

Fig. [T] shows the adiabatic eigenfrequency differences 
between all the possible aforementioned options, only ra- 
dial modes are depicted. T he equilibrium mo del used is the 
last of the stepl of Task 2 CMoya et alil2007h (this volume), 
with 4042 mesh points, 1.5Mq and Xc = 0.4. As reference 
we have calculated eigenfrequencies with the following op- 
tions: X ={£ = 0, no Richardson, p', dp = 0, Inr). For each 
comparison we have changed one single degree of freedom, 
remaining the rest unchanged. We show the differences ob- 
tained in the range [200,2500] /iHz, that is, from the fun- 
damental radial mode to a frequency slightly larger than the 
cutoff frequency (around 2250 /iHz). Top panel presents the 
differences obtained when two outer mechanical boundary 
conditions are used. Reference line is dP = 0, and the com- 
parison is with the use of the isothermal reflective wave outer 
boundary condition. The differences for large frequencies 
are of the order of units of /iHz. These differences are sim- 
ilar to those obtained by other codes (J.C. Suarez, private 
communication), but a larger study of these differences is 
still needed. The use of higher order integration procedures, 
as the Richardson extrapolation, do not change significantly 
this differences. In bottom panel, the rest of the compar- 
isons are depicted. The Richardson extrapolation gives dif- 
ferences of the order of tenths of /iHz, the use of r/P as 
integration variable provides small differences always lower 
than 0.008/iHz. The LAWE differential equations and the 
use of the Lagrangian variation of the pressure as eigenfunc- 
tion provide similar differences smaller than 0.05/iHz, but 
its profile is no t constant. For a com prehensive study of these 
differences see lMova et alj (120071) (this volume). 



3 Non-adiabatic resolution 

Additional information for asteroseismology is provided by 
GraCo with the resolution of the non-adiabatic set of dif- 
ferential equations described in Unno et al. ( 1989) (pp 261 
ff). In the non-adiabatic resolution the eigenfrequencies and 
the eigenfunctions are no longer real. This makes it possi- 
ble to obtain the so called non-adiabatic observables: 1) The 
"phase-lag" {^^ = 0(^;.) — (/)(r)) defined as the phase differ- 
ence between the Lagrangian variation of the effective tem- 
perature and the radial displacement. 2) The relative varia- 
tion of the Lagrangian variation of the effective temperature 

( T-*"-^-^ ). And 3) The energy balance of each mode measured 
by the "growth rate", directly related with the imaginary part 
of the modal eigenfrequency. 
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Fig. 1 Frequency comparison of the different adiabatic resolutions 
given by GraCo (see text). Only radial modes are shown. The equi- 
librium model used is a l.SMg and Xc = 0.4 model, with 4042 mesh 
points. 



The code uses the adiabatic solutions obtained for a given 
mode as trial functions for the non-adiabatic relaxation pro- 
cedure. The i nner b oundary conditions are those described in 
lUnno et alTI ( Il989t) (p. 229). It is in the outer region where 
GraCo presents some complexity. The code can treat the 
photosphere as a boundary condition or introduce the atmosphere- 
pulsation interaction resolution described by 
iDupret et al. I (l2002h . This interaction is described imposing 
the atmosphere to be in thermal equilibrium and the diffu- 
sion approximation for the radiative flux to be no longer 
valid. In this case two different sets of differential equations 
are solved, one for the stellar core and envelope and another 
for the atmosphere. A transition layer and outer boundary 
conditions for the atmosphere must be defined. 

Fig. |2] shows the values of the non-adiabatic observables 
for a standard 5 Scuti model of l.SMp,. In this figure we 
can see how the values in the case labeled "with" (where the 
atmosphere-pulsation interaction is included) are clearly dif- 
ferent from those labeled "without" (the photosphere treated 
as the outer boundary layer). This illustrates the importance 
of the inclusion of the atmosphere -pulsation interaction for 



Granada oscillation Code (GraCo) 













vio,t^ 












16 






























' 






' 


' 






X 
















with + 




14 


'\ 
















without X 


- 


12 




















- 


h" 


- 




















\ 10 

CO 


* \ 




















8 
6 


*^ 


'"x 


^^Xx 


^ 




; 


X 


XXX 


XX ^ XX 


_l. 



0.01 



0.02 0.025 

Q (in days) 

l.S Mq, a=l, CEFF, X^=0.44 





>>< 




' 






' ' 






>>( 














190 




>!- 








with + 
without X 




180 






X 










I/O 


+^ 




'x 










160 




+^ 


x^ 








" 


ISO 




+ 




^ 






- 


140 






+ 




X 




,; 


130 






+ 








- 


120 

1 1A 






++ 


- .., , 


+ 


+ ++ ^ 


" 



0.01 0.015 0.02 0.025 0.03 0.035 

Q {in days) 

Fig. 2 |5Teff/Teff| (top panel) and (^'^ (in degrees, bottom panel), as a 
function of the pulsation constant Q (in days) for different modes with 
spherical degrees / = 0, 1,2,3. A model of a l.SMp, is studied with 
Xc = 0.44, a MLT parameter a = I and the CEFF equation of state. 
Results obtained "with" (+) and "without" atmosphere (x) in the non- 
adiabatic treatment are compared. 



a better description of the non-adiabatic observables. When 
the atmosphere is treated as a boundary condition, not all the 
heat exchanges here are correctly taken into account, there- 
fore the phase-lags obtained are closer to the adiabatic pre- 
diction (180°). On the other hand, the system of equations 
modeling the atmosphere-pulsation interaction takes into ac- 
count these non-adiabatic processes in the atmosphere, and 
the resulting phase-lags are much smaller than 180°. This 
atmosphere-pulsation interaction has not influence upon the 
growth rate, since it takes place in layers not relevant for the 
driving of the modes due to their very small density. 

All of the above-mentioned calculations have a direct in- 
fluence on the phase difference - amplitude ratio diagrams 
used to discriminate oscillation modes, phase-lags, as well 
as relative variations in |5Teff/Teff| and dge/ge (also cal- 
culated in the non-adiabatic resolution) can be used to over- 
co me the uncerta i nties i n previous phase-ratio color diagrams. 
In lGarrido et al. I (|199C^ these discrimination diagrams were 
made using parametrized values for departures from adia- 
baticity and phase lags. The only remaining degree of free- 



dom is now the choice of the MLT a parameter in order to 
describe the convection. Therefore, discrimination diagrams 
depend only on this parameter, as is shown in Fig. [3] for 
the same equilibrium model used in fig. |2] Theoretical pre- 
dictions are plotted for two specific Stromgren photometric 
bands {{b — y) and y) using three different MLT a parame- 
ters in the fundamental radial mode regime (pulsation con- 
stant near 0.033 days) and in the 3rd overtone regime (near 
0.017 days). 

A clear separation between the /-values exists for periods 
around the fundamental radial. Similar behaviour is found 
for other modes in the proximity of the 3rd radial overtone, 
although for these shorter periods some overlapping start to 
appear at the lowest /-values. They also show the same trend 
as for the fundamental radial mode: high amplitude ratios 
for low MLT a and the spherical harmonic 1 = 3. 



4 Conclusions 

GraCo is a complete software for the resolution of systems 
of differential equations related with the stellar oscillations. 
The general integration scheme used in t he code is the Henyey 
relaxation method as it is explained in lUnno et al. I (|T989). 
In the adiabatic frame, different integration schemes can be 
used to obtain the eigenfrequencies and eigenfunctions: 1) 
Two outer mechanical boundary conditions, 2) two choices 
for the set of eigenfunctions, 3) two choices of the integra- 
tion variables, 4) the use or not of the Richardson extrapo- 
lation and, 5) for radial modes the use of the LAWE second 
order differential equation or set £ = in the standard non- 
radial differential equations. 

On the other hand, the main characteristic of the code is 
that the non-adiabatic set of differential equations can also 
be solved. Two different treatments of the photosphere can 
be used. The first considers the atmosphere as a single bound- 
ary layer, and the second describes the atmosphere -pulsation 
interaction. This makes it possible to obtain more accurate 

non-adiabatic observables (phase-lag and -j^) crucial for 

modal identification through the multicolor photometry. On 
the other hand, non-adiabatic studies also allow to study the 
modal energy balance, giving a theoretical range of over- 
stable modes. 

Finally we want to remark that the analytical expressions 
of the differential equations, eigenfunctions and boundary 
conditions, and the numerical proc edure followed in GraCo 
can be found in lUnno et al. I (Il989h . 
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